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DISCUSSION OF A DIRECT STEP METHOD FOR 
COMPUTING WATER-SURFACE PROFILES 
PROCEEDINGS-SEPARATE 180 


I. H. Sremvserc.t—The advantages of any graphical or semigraphical 

4Superv. Hydr. Engr., Corps of Engrs., U. 8. Dept. of the Army, San Francisco, Calif. 
solution over the trial-and-error methods of flow line computations lie in the 
frequent use of the derived curves. Where only a relatively few profiles of 
given discharges and conditions are required, it has been found practicable to 
draw the elevation-area curves and elevation-(A R!) curves to a convenient 
scale and use the trial-and-error method for tracing the water-surface profile. 
There have been many occasions, however, when it was necessary for the 
writer to compute numerous water-surface profiles over the same stretch of 
river. It would have been an insurmountable task were it not for the fact 
that it was possible to use graphical solutions. In most of the graphical 
methods corrections for velocity head were not considered. 

The author has presented a relatively simple graphical method of computing 
directly the water-surface profiles in natural and artificial watercourses. 
Where profiles for discharges of different magnitudes are required, it would be 
desirable to present the pairs of curves for the discharges on separate axes for 
each reach. Thus, only the curves of Z + $(Z) for station 80 + 00 and the 
curves of Z + f(Z) for station 85 + 00 for the discharges of 33,500 cu ft per 
sec and 40,000 cu ft per sec would be shown in Fig. 2. The curves for the 
other reaches could be shown in a similar manner. 

Mr. Ezra uses a weighted value for the velocity of flow when both channel 
flow and overbank flow are involved. A suggested method is to weight the 
velocities in accordance with the respective channel and overbank discharges. 
A weighted-area curve which can be applied directly to the discharges to obtain 
the weighted velocities can then be derived. Therefore, 


Ve +Q Vo 
v 


in which P is the weighted velocity in feet per second, Q. is the channel flow 
in cubic feet per second, V. denotes the velocity of flow in the channel, 
Q. represents the overbank flow in cubic feet per second, V, is the velocity of 
overbank flow, and Q is the total discharge in cubic feet per second. 

Since 


then 


in which A, is the weighted cross-sectional area and A, and A, are the cross- 
sectional areas of the channel flow and the overbank flow, respectively. 
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shy 


From the Manning formula, 


and 
Q = (C. + C.) hi, = Ch, 
in which 


1.486 A 
n 


C= 


and h, is the friction loss in feet in the distance L. From Eq. 68, 


C? h?; hy + hy; 
Aw A. 


from which 


Either Eq. 9b or Eq. 9c can be used to compute the weighted-area curve. 
Where the channel and overbank lengths are equal, the term L will cancel and 
therefore it can be omitted from Eq. 8. Where the reach lengths are different 
for channel and overbank, the average lengths for adjacent reaches upstream 
and downstream from the section should be used in Eq. 7. F 

If the velocity head rather than the velocity is weighted, the expression 
for the weighted-area curve will be 


1 
Cc Cc A, 
A number of graphical and semigraphical charts and graphs have been 
developed for eliminating the trial-and-error computations in tracing the 
water-surface profile. The writer devised a method which required a nomo- 


graph for the solution.’ In most methods the velocity-head corrections have 


5“*The Nomograph As an Aid in Computing Backwater Curves,” by I. H. Steinberg, Civil Engineering, 
Vol. 9, 1939, p. 365. 


not been taken into account. However, Francis F. Escoffier has presented a 
solution which is similar to the direct step method developed by Mr. Ezra 
in the respect that pairs of curves are plotted for each reach. A sloping 
straight line drawn from the known point on one curve to the intersection of 
the second curve gives the required water-surface elevation. The slope of the 
line is a function of the discharge. Thus, only one set of basic curves is 
required for a wide range of discharges. 
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AS 
or 
Ay = A. 
1+ (4:)() 


The writer has also developed a graphical solution utilizing “slope curves” 
in which velocity head corrections are provided for. In this method the 
water-surface elevations at the lower end of the reach and the discharges are 
plotted as ordinates and abscissas, respectively, whereas the falls, or differences 
in water-surface elevation between the upstream and downstream end of the 
reach, are plotted as the parameter. The curves can also be plotted with the 
elevation at the upper end of the reach as the parameter. Eq. 1 can be written 
as 

v3, 
Zi+ Z2+ 29 
in which 


V%, 


Using the Manning formula for computing the friction loss, 


in which 


1486 An R},, 
n 


The subscript m is used to denote the mean value for the reach, and atten- 
tion is directed to the different method of applying the hydraulic elements to 
obtain the average slope in the reach from that used in Mr. Ezra’s derivation. 
The mean value for the reach can be based on the average of the end value 
or on the value as defined by the midpoint elevation of the reach. The writer 
uses the latter system because of the greater simplicity in computing the 
basic slope curves. 

Eq. 1 can also be written as 


However, 


Therefore, 


2 


= = Qn on _ Ps) 
Zi (2) 29 +K 29 29 ... (18a) 
and 
Q 
and 
r=(2) - 4,) |“ K)} = (4) 
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Water -Surface Elevation, in Feet 


Wate -Surface Elevation, in Feet 


1 
18 20 22 24 26 28 30 8 30 32 34 36 38 40 42 44 46 
Value of ARix ee Area, in Square Feet x 107? 


Fra. 4.—Errvation—(A R!) Crrves Fic. Crrvrs 


in which 


Any given values of midpoint elevation and of fall, F, in a reach will define 
corresponding values of C,,, Ai, and A». The value of M can be computed 
from Eq. 15a and the value of Q then determined from Eq. 156. For develop- 
ment of the “slope curves,”’ Eqs. 15a and 156 are solved for a range of values of 
midpoint elevation and fall. The resulting points are plotted as curves on 
appropriate graphs. 

To illustrate the application of the “slope curve” graphical method in 
which a velocity-head correction is included, Example No. 1 will be used. 
Values of K equal to 0.5 for decelerating flow and K equal to — 0.1 for ac- 
celerating flow have been adopted. The negative value of K for accelerating 
flow has been used to conform with Eq. 15a. 

Table 7 contains the same basic data as given in Table 1, except that the 
values of A R! have been computed. Elevation-(A R!) curves and elevation- 


Water- Surface Elevation 
at Station 80 +00, in Feet 
Water-Surface Elevation 
at Station 85 +00, in Feet 


Discharge, in Cubic Feet per Second x 10~? Discharge, in Cubic Feet per Second x 107? 


Fia. 6.—Stopr Crrves For Rracn A Fic. 7.—Stopr Crrves ror Rracn B 


23 = 
21 21 
Lf 
18 18 Pa | a 
| 
(1— K)........(15a) | 
and | 
22 22 
|| 
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Water -Surface Elevation 
at Station 89 + 00, in Feet 


~ 


25 30 35 40 45 50 
Discharge, in Cubic Feet per Second x io 


Fic. 8.—Stopr Crraves ror Rracn C 


area curves, based on Table 7, are shown in Figs. 4 and 5. Table 8 illustrates 
the computations for reach A (station 80 + 00 to station 85 + 00). The 
values of A R! (Col. 5) and areas (Cols. 6 and 7) were obtained from the 
curves in Figs. 4 and 5. The resulting slope curves are shown in Fig. 6. The 
curves for reaches B (station 85 + 00 to station 89 + 00) and C (station 
89 + 00 to station 94 + 00) are shown in Figs. 7 and 8. 


TABLE 7.—Hypravtic Properties oF Cross SECTIONS 


Water-surface Cross-sectional Wetted perimeter 


Station elevation, Z, area, A, in ARI 
in feet square feet P, in feet 
(2) (3) 


(4) 
238 
248 
258 


| 
| 
YF. Wa W 
| -9 o 
| 
18.0 3,000 16,210 
80 + 00 20.0 3,500 20,400 
22.0 3,970 24,580 
: 18.0 3,580 265 20,400 
85 +00 20.0 4,100 273 25,000 
22.0 4,630 282 29,800 
18.0 3,080 248 16,400 
89 + 00 { 20.0 3,550 258 20,500 
22.0 4,050 268 24,750 
18.0 2.650 225 13,710 
94° + 00 20.0 3,100 235 17,350 
22.0 3,550 245 21,000 
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TABLE 8.—CompuTaTIONS FOR DEVELOPING THE 
Store Curves For Reacu A 


Warer-Surrace Cross-SecrIONAL 
ELEVATION, IN AREA, IN SQUARE 
in feet | in feet Zs. Qs 
Station Station Station Station y 
80 +00 | 85 + 00 80 +00 | 85 +00 


(3) (4) (S)« (6) (9)4 (10)¢ 


27,800 
30,000 
31,800 
33,600 
35,050 
36.400 


Discharge, 


31,050 
33,500 
35,600 
37,400 
38,900 
40,500 


wine 
BS | 


| 
oN 
| 


37,050 
39,950 
42,400 
44.500 
46,600 
48,200 


0 
.2 
4 
6 
8 
0 


@ Values of 4 R} obtained from Fig. 4.» Values of Cm determined from Eq. 126 with n = 0.030 and 
L = 500 ft. © Values of A obtained from Fig. 5. 4 Values of M? computed from Eq. 15a with K = 0.5 
when Ai < 4: and K = —0.1 when Ai > Ax. © Values of Q determined from Eq. 150. 


To compute the water-surface profile, it is necessary only (1) to enter the 
slope curves with the given elevation at the lower end of the reach and the 
given discharge, (2) to read the corresponding fall, and (3) to add to the given 
elevation to obtain the elevation at the upper end of the reach, which then 
hecomes the given elevation at the lower end of the next upstream reach. This 
procedure is illustrated in Table 9. The computed water-surface profile 


TABLE 9.—ComputTaTION OF WATER-SURFACE PROFILE FOR A VALUE 
oF Q = 33,500 Cu Fr per Sec 


Water-surface elevation, 
Station Reach Z, in {cet Fall, F, in feet 


(1) (3) 


80 + 00 18.50 
85 +00 19.83 
89 +00 20.04 
94 +00 20.52 


shown in Table 9 agrees closely with that shown in Table 6. 
It should be noted that the slope curves have been drawn for computing 
upstream. The same curves can be readily adapted, however, for computing 
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| 
{ 
(1) (2) | 
‘ 1.0 18.0 19.0 3,000 2.39 
12 17.9 19.1 2'970 247 
14 178 19.2 255 
18.5 16 17.7 193 19,450 | 43,100 | 4 9'999 2.63 
18 176 194 2'900 271 
20 17.5 195 2'880) 2.80 
1.0 19.0 20.0 3,255 2.37 
12 18.9 20.1 3'230 2.45 
14 18.8 20 2 3/205 254 
19.5 16 18.7 20.3 21,600 | 47,800 | 4 3'180 2 62 
18 18.6 20.4 3.155 271 
20 18.5 20.5 3.130 279 
10 | 200 21.0 3,500 2.40 | 33,900 ) 
12 19.9 21.1 3.475 248 | 36.600 
14 19.8 21.2 3.450 256 | 38.900 
20.5 16 19.7 21.3 ¢ ° | 23-800 | 52,600 | 4 3°495 265 | 40.900 
18 19.6 214 3400 | 4470 | 274 | 42,900 
20 19.5 215 3.380 | 4500 | 282 | 44.400 } 
1 21.0 22.0 3,740 | 4.630 | 241 
I 20.9 221 3715 | 4655 | 249 
1 20. 22.2 3.690 | 4.685 | 258 
21.5 1 20.7 223 26,000 | 57,500 | 4 3'665 | 4.710 | 267 
1 20.6 22.4 3.640 | 4.735 | 275 
2 20.5 22.5 3.620 | 4.760 | 283 
| 
| 
(4) 
1.33 
0.21 
0.48 


downstream by trial and error. In this ease—(a) the elevation at the lower end 
of the reach is assumed and (6) the fall is obtained from the slope and added 
to the assumed elevation. This elevation is checked against the given elevation 
at the upper end of the reach. Usually two, or perhaps three, trial computa- 
tions are reyuired to obtain the correct fall. An illustration of the procedure 
is given in Table 10. 


TABLE 10.—DownstreamM ComputaTION OF WaTER-SURFACE PROFILE 
For A Vaturt or Q = 33,500 Cu Fr per Sec 


Computed 
wafer-surface 
elevation at 
upper end of 
reach 


Trial water- 
surface elevation Fall, F, 
at lower end in fect 
of reach 


Water-surface 
Station elevation, Z, 
in feet 


(4) (5) (6) 


20.60 0.42 21.02 
20.57 0.43 21.00 


20.35 0.19 20.54 
20.38 0.19 20.57 


19.25 1.18 20.43 
19.20 1.18 20.38 


C. A. M. Gray,® A. M. ASCE.—An ingenious variant of the standard step 

® Senior Lecturer in Civ. Eng., Univ. of Sydney, Sydney, Australia. 
method has been presented by the author. As the solution is direct and the 
labor involved is considerable, the question of how long the steps are to be made 
is most important. It seems worthwhile, therefore, to investigate the measure 
of the error introduced by the selection of an arbitrary length of step. 

Eq. 2, neglecting the eddy-loss term, is obtained from the ‘“‘exact’’ integral 


relationship, 
V3, f 
22+ =Z 


when the integral is evaluated by the trapezoidal rule. An alternative value 
would be Sy L in which Sy, is the value of S at the midpoint of the interval L. 
If the integral is denoted by /, it can be shown that 


>1> Sub 


Sy L>I>—>— 


depending on whether the curve of S versus z is concave upward or concave 
downward. In the former case, Mr. Ezra’s solution would give too small a value 
for Z, and in the latter case it would result in too great a value. Using the mean 
value of S will give an opposite result. An estimate of Sy can be obtained by 
use of the data required for Fig. 2. If 


y2 
(z 


(1) (2) (3) 
89 +00 20.57 
B 
8 +00 20.38 
Cc 
80 + 00 19.20 
or 
V3, 
and 
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in the range 0 < x < L, then 


Hence, 
29 2 


V3, AL 


+ 


in which 


By = (+2) 
29 zeL/2 


From the properties of the section at z = L/2, limits ean be set for Sy. Un- 
fortunately, the data are not available. For purposes of investigation it can 
be assumed, however, that the cross section is constant between stations. 
From Tables | and 3 it can be seen that, in the interval from station 80 + 
00 to station 85 + 00, 
0.00076 < S < 0.00162 


and the value of E at station 80 + 00 is 20.31. Hence, from Eq. 19a, 
0.40 + 20.31 = 20.71 > Ey > 0.19 + 20.31 = 20.50 


From Table 1 
0.00149 > Sy > 0.00139 


S, S, 


which is greater than the mean value 2 


, indicating that the curve of S 


versus x is concave downward. 
Hence, at station 85 + 00, 


E < 20.31 + 0.75 = 21.06 
and from Table 1 
Z < 20.02 


Under these conditions, at station 85 + 00, 


19.83 < Z < 20.02 
and 
20.90 < E < 21.06 


Upper and lower bounds to Z at station 85 + 00 have thus been established, 
and the separation of these values gives a measure of the accuracy achieved 
with the assumed length of step. 

Bos Buexwer,’ A. M. ASCE.—A method has been presented by Mr. Ezra 

7 Hydr. Engr., TVA, Knoxville, Tenn. ; 
which should be welcomed by those engaged in computing water-surface 
profiles in natural streams or reservoirs. If profiles for several flows and down- 
stream starting elevations are to be determined, use of this method, as opposed 
to the standard trial-and-error procedure, will save time and avoid drudgery. 
The method has been tried successfully on a typical stream, and it is believed to 
be superior to most others. 
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This discussion (1) seeks confirmation or simplification of a method pre- 
sented by the writer for the consideration of overbank flow ; (2) offers alternate 
plotting procedures; and (3) seeks further confirmation of (a) the application 
of channel-roughness coefficients which vary from reach to reach, (b) the de- 
termination of these coefficients from observed flood profiles, and (c) the 
determination of channel coefficients under conditions in which there is a 
substantial amount of overbank flow. 

Inclusion of Overbank Flow.—Tabie 11, showing hydraulic properties, is a 


TABLE 11.—Hypravuuic Properties oF Cross Secrions ror Q = 33,500 
Cu Fr per Src, n for CHANNEL = 0.03, AND n FOR OVERBANK = 0.05 


Friction gradient 


Water-surface 
elevation Z, 
Cross-sectional 
area .1, in 
Hydraulic radius 
R, in feet 


in feet 
Weighted V2 


square feet 
Velocity V, 


(3) 


=~ 


(10) qt) 


0.0416} 0.00173 


1,012,000 2,800,000 
30,900 92 3,890 


0.0321) 0.00103 2,803,890 
1,891,000 
15,620 


1,315,100 0.0255} 0.000650 1,906,620) 
1,005,000 
00 0 


0.0333) 0.00111 
1,237,000 8.09 | 2,170,000 
12,660 339 433 


1,249,660 A 0.0268) 0.000718 2,170,433 
1,482,000 d 1,502,000 
78,100 5,800 


1,560,100} 33,500 | 0.0215] 0.000462 1,507,800 


variation of Table 1 in that flow in an assumed overbank has been added. 
For brevity, only the reach from station 80 + 00 to station 85 + 00 has been 
used. Table 11 shows, for each elevation, two values for the area and two 
values for the hydraulic radius. The upper figure is the channel value used in 
Table 1, and the lower tigure is the overbank value that has been assumed for 
illustration. The value of n for the overbank has been taken as 0.05. If 
discharge is represented by 


1.486.A Ris! 
n 


1.486 A RI 
n 
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Si 


. . (236) 


N 
§ s 
= 
(1) (2) |_| (4) (S) (6) (7) (9) ee (12) 
‘3,000 | 12.60] 305,000 
18.0 { 0 | 0.00 0 0 
33,500 | 11.17 1.94 
(3,500 | 14.11 
80+00 | 4200 | | 500 | 3.00 
| 1,042,900] 33,500 83.7 1.30 
22.0 ‘panes 6.00 98,100} 2,500 
| 6.9 0.88 
18.0 0; 0 
85 +00 20.0 300 1.70 
64.8 1.01 
22.0 900 | 5.00 
4° () 0.70 
then 


In Col. 5, Table 11, Q S~! can be computed first for the channel section, and 
next for the overbank; then the two values are added. For El. 20.0 at station 
80 + 00, Q S! is 1,012,000 in the channel, 30,900 in the overbank, thus total- 
ling 1,042,900. The true total discharge is known to be 33,500 cu ft per sec. 
Therefore, Cols. 7 and 8 (Table 11) ean be determined. With S! known, the 
parts of the total flow carried respectively by the channel and overbank can be 
computed and entered in Col. 6, Table 11. The velocities and weighted velocity 
head for the total section can be found and entered .n Cols. 9 through 12, 
Table 11. 

Table 12 shows, for this same reach, the computation for @(Z) and f(Z) for 


TABLE 12.—Compuration ror Z + f(Z) anp Z + $(Z) ror Reacu 
FROM Sration 80 + 00 to Srarion 85 + 00 wirH 7 FoR 
CHANNEL = 0.03, n For OveERBANK = 0.05, 

AND L/2 = 250 Fr 


Vatves or @(Z) or f(Z) | Vatves or Z + on 
ror Q, In Cr Z+(Z) vor Fiow Q, 
Fr per Sec in Cu Fr per Sec 


Vatves For Q = 33,500 
Cu Fr eer Sec 


Weighted 
10,000) 20,000) 30,000) 40,000 | 10,000 | 20,000) 30,000 | 40,000 
0.43 1.94 2.37 0.21 | 0.84 § 3.37 | 18.21 | 18.81 | 19.90 | 21.37 
0.26 1.30 1.56 0.14 a 23 2.22 | 20.14 | 20.56 | 21.25 | 22.22 
0.16 0.88 1.04 0.09 4 . 1.48 | 22.09 | 22.37 | 22.83 | 23.48 


0.28 1.36 0.10 : } 1.54 | 18.10 | 18.38 | 18.87 | 19.54 
0.18 1.01 0.07 1.18 | 20.07 | 20.30 | 20.66 | 21.18 
0.12 0.70 0.05 0.83 | 22.05 | 22.21 | 22.46 | 22.83 


Water-surface 
elevation Z, 


in ft 


Station 


the base flow of 33,500 cu ft per see using the slopes and weighted velocity heads 
found in Table 11. Other entries in Table 12 will be explained subsequently. 

Alternate Plotting Procedures.—Adaptation to specific needs as well as to 
individual preference determines the most convenient method of plotting. For 
example, if curves were to be prepared which would be suited to a wide variety 
of flows, the form used by the author would become unwieldy. It would seem 
logical, therefore, to use one plotting sheet for each reach, showing on each a 
family of discharge curves. An example for the reach from station 80 + 00 
to station 85 + 00 is shown in Fig. 9 which was prepared from Table 12. 
Values of @(Z) and f(Z) in Table 12 for the various flows were computed from 
the base value for 33,500 cu ft per sec, using the ratio of the discharges squared, 
as suggested by Mr. Ezra. The solid lines in Fig. 9 were drawn from computed 
values of Z + $(Z) and Z + f(Z). The dashed intermediate lines were inter- 
polated by use of construction lines based on the computed curves. The form 
of this plot differs from that presented by the author. In this form, if eddy 
losses are ignored, the upper curves can be entered with the known elevation 
at the downstream section and flow—for example, Fl. 20.32 at station 80 + 00 
and 37,000 cu ft per sec. This intersection is carried vertically downward to 
the lower curves until the proper discharge is reached. This lower intersection 
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f 
{ 
80+00 | {20.0 
\22.0 
(18.0 
85+00 | (20.0 
(22/0 
f 


gives a resulting elevation of 21.28 at station 85 + 00. The values of Z + $(Z) 
and Z + f(Z) need not be tabulated since they only represent guide lines 
connecting the upper curves for station 80 + 00 to the lower curves for sta- 
tion 85 + 00. 


The basic form of the curve in which Z is plotted against Z + $(Z) and 
Z + f(Z) has serious limitations if the water-surface elevations at a given sec- 
tion will vary over a wide range. A scale on a convenient size of graph paper 
can be read with sufficient accuracy if the graph covers an elevation range of 
perhaps 10 ft, but not if it covers a range of 50 ft or more. For any given 
reach, however, there usually is a reasonable limit in the elevation differences 
at the two ends. The functions ¢(Z) and f(Z) are closely related to the rise 
in a reach, suggesting a rearrangement of Eq. 4. Neglecting eddy losses, 
Eq. 4 can be written as 


Z, = Z2 + $(Z2) — f(Z1) 


For any reach, a curve can be drawn for the downstream section showing 
@(Z) versus Z. Another curve can also be drawn for the upstream section 
showing f(Z) against Z+ f(Z). Fig. 10 is an example of this in which the 
reach from station 80 + 00 to 85 + 00 has been used. Fig. 10 was plotted 
from the data in Table 12. Table 13 shows the water-surface profile for a flow 


TABLE 13.—Warer-SurFace 
PROFILE FOR Q = 37,000 
Cu Fr per Sec 


Water- 
surface Value | Value | Value 


Station elevation of 
o(Z) 
in feet 

(b (2) (3) (4) (5) 


80 + 00 20.32 1.78 | 22.10 
85 + 00 21.28 


of 37,000 cu ft per see and illustrates the use of the curves. Starting with 
El. 20.32 or Z at station 85 + 00, ¢(Z) (Col. 3, Table 12) is found to be 1.78 
from the upper curves. This value is added to Z to get a value of Z + $(Z) 
equal to 22.10 (Col. 4, Table 13). From Eq. 4, however, this value equals 
Z+f(Z) at station 85 + 00. Entering with this value of 22.10 in the lower 
curves, f(Z;) (Col. 5, Table 13). is found to be 0.82. The value of Z at station 
85 + 00 (El. 21.28) is Col. 5, Table 13, subtracted from Col. 4, Table 13. In 
this procedure, elevations are not read from the curves. Instead, factors com- 
prising the difference in elevations are read. The operation requires more 
tabulating than when Fig. 9 is used; but with the same size graph paper, it is 
far more accurate whenever wide elevation ranges must be accommodated. 
Of course, it is presumed that reaches will be short enough so that values of 
¢(Z) and f(Z) will be sufficiently small to plot with reasonable accuracy. 
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Elevation at Station 80 +00, in Feet 


WN 
\ 


§ 
+ 
wo 
c 
2 
= 
2 
= 
2 
= 
> 
2 
w 


Fro. 9.—Crrves For WaTeR-ScRFACE PROFILE 


Roughness Coefficients and their Determination —Whenever possible, it has 
been the writer’s practice to use roughness coefficients that were computed 
from observed flood profiles for known flows. If a choice of floods for this 
purpose is available, the preferred flood is one which just fills the channel part 
of the cross section without extending into the overbank. The resulting com- 
puted channel coefficients are then tested on one or more higher floods to aid 
in assigning a suitable roughness coefficient to the overbank. Channel- 
roughness values usually vary enough, and are critical enough, so that different 
values are required in each reach. Overbank values are (a) more difficult to 
compute with assurance, (b) are usually less critical than channel values, and 
(c) are often assumed to be constant for long stretches of a stream. Channel- 
roughness factors are computed, and applied, with the friction slope based on 
the average of the areas and the hydraulic radii at the two sections bounding 
each reach. 

Mr. Ezra uses a different procedure in the application of area and of hydraulic 
radius. The rate of friction slope is found for each section using its area and 
hydraulic radius. The total friction slope in a reach is, then, the average of the 
rates at its two end points times the length. If roughness factors vary in each 
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reach, the rate of friction loss at each section must be computed twice—once 
for use in an upstream direction, and once for use in a downstream direction. 
This would require two columns for S in Table 1, one headed S, and one Sq. 
It would be similar to the need for both Cols. 3 and 4, in Table 2, in which } L 
may differ in the upstream and downstream directions. Also, Col. 3, Table 2, 
would then be entitled } Sz La and Col. 4, Table 2, would be entitled 3 S, Ly. 

When overbank flow is not involved, roughness factors can presumably be 
computed from a marked flood profile of known flow using Eq. 3. All quantities 
in this equation are known except the value of n. When n must be computed 
from a flood involving appreciable overbank flow, however, the determination 
of the channel values, even with the value of n in the overbanks assumed, 
becomes a cumbersome operation. Were it not necessary to weight the channel 
and overbank velocity heads, the overbank flow could be computed first and 
deducted from the known total flow, thereby leaving the necessary known 
quantities with which the channel roughness can be computed. This would 
be a reasonable expedient if the overbank flows were not large. When overflow 
is large, however, and when it appears necessary to use a weighted velocity 
head (even with an assumed overbank n), it seems necessary to make successive 
estimates of the value of n for the channel until a value is found that will satisfy 
Eq. 3 (modified to include overbank flow). An application of this principle 
to a typical stream required the use of a 20-eulumn table to reach a satisfactory 
solution. 

Monrr M. Kanson.*—The method proposed by the author is a valuable 

* Asst. Prof., Faculty of Eng., Univ. of Alexandria, Alexandria, Egypt. 


contribution to the problem of plotting water-surface profiles. It facilitates 
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plotting and saves time, especially once the necessary tables and graphs are 
computed and prepared. 

Mr. Ezra describes his method as being a general one. Actually, however, 
the method is limited to eases involving gradually varying flow in which de- 
parture from hydrostatie conditions is so negligible that the potential energy 
can be expressed by the terms Z, or Zs in Eq. 1. The method is based on the 
assumption that the velocity is uniformly distributed over the area of any cross 
section. 

The term “S”’ varies with the rugosity factor which can be supposed to be 
constant betweea two consecutive stations, as well as with the velocity of flow 
and with the proportions of the section of the channel. The average rate of 
friction loss of head can then be taken to he equal to 


g - ( _Vavwe 
1.486 Raver®!? 


Assuming, in a uniform eanal, that 


2 R, + 


n? (V, + 7.)? 


(3) 2.2 (R, ( 7) 
Eq. 27 shows that Sayer, as derived, does not equal } (S, + S2), nor 
n? V2, V2, b h 

} 32 \ Rus, + Ri, as proposed by the author. 

It should be realized that the area, the perimeter, the hydraulic mean depth 
2 
or radius, and consequently a and ae are functions of the depth of 
flow and not of the elevation of the water surface above a chosen datum. It is 
not correct, then, to denote that (in general) 


2 L 2 n2 
2 Ra = (2) and f (2). 
This is true only in particular cases for particular cross sections selected before- 
hand. 

It would be helpful if a plot were presented of the same surface profiles 
computed by Mr. Ezra’s method and according to the methods developed by 
Bresse, Bakhmeteff, and Tolkmitt (as well as some of the simplified analytical 
and graphical methods) to show the discrepancy. 

Arruour A. Ezra,’ J. M. ASCE.—Mr. Kansoh is correct in stating that the 

* Senior Harbor Engr., Port of San Francisco, Calif. 
proposed method is limited to conditions in which the potential energy can be 
expressed by the terms Z; or Z. in Eq. 1. This limitation corresponds to the 
boundary conditions (found in the theory of varied flow in open channels) on 
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which are based the generally used methods of determining water-surface 
profiles. The velocities used are usually the mean or weighted velocities at a 
section. The area, the perimeter, the hydraulic radius, and consequently 
the terms a and soa are functions of the depth of flow and the eleva- 
tion of the water surface above a chosen datum. The reason for this inter- 
relationship is that the depth of flow is a function of the difference between the 
water-surface elevation and the elevation of the channel bottom. The eleva- 
tion of the channel bottom is expressed implicitly in the values of A and R that 
correspond to a given water-surface elevation at a section. 

Mr. Kansoh’s question regarding the value of S,ver appears to be a matter 
of definition. The writer defined S,,., as that value of S which is the arith- 
metic mean of the values S,; and S». These values of S are determined indi- 
vidually from Manning’s formula using V;, Ri, and V2, Re, respectively. Mr. 
Kansoh defines S,,y.r as that value of S which is obtained by substituting in 
Manning’s formula the arithmetic mean of V; and V,. together with the arith- 
metic mean of R, and Ry. Naturally, the value of Sayer as defined by Mr. 
Kansoh will not equal } (S; + S:). The writer’s usage of 3 (S, + S,) as the 
average rate of loss of head corresponds to that used by George E. Russell.” 

1” “Hydraulies,"’ by George E. Russell, Hemy Holt and Co., New York, N. Y., 5th Ed., 1947, p. 292. 
The writer’s method gives plots of surface profiles corresponding to those ob- 
tained by the standard step method which is based on trial and error. Mr. 
Kansoh’s suggestion of plotting the same surface profiles for comparison by 
Bresse’s, Bakhmeteff’s, and Tolkmitt’s methods is profitable. 

The writer is indebted to Mr. Gray for his analysis of the relative accuracy 
of the approximation used to establish the friction-head loss in a reach in 


comparison to the results obtained by using f Sdzxand Sy L. From Table 


3 the difference in water-surface elevations between station 85 + 00 and 
station 80 + 00 is 1.33 ft. Using Mr. Gray’s method of analysis, the difference 


2 
in water-surface elevations computed by using f S dz is between 1.33 ft and 
1 


1.52 ft. On the basis of this analysis, the writer’s approximation yields a result 
for a 500-ft-long interval that is within 12% of the result that would be obtained 


by using fs dx. This inaccuracy is small in comparison to the inaccuracies 


involved in the choice of a value of the roughness coefficient n, especially since 
S varies as n?. 

Mr. Buehler has made a valuable contribution to the method of computing 
water-surface profiles. His suggestions for alternate plotting procedures are 
very appropriate. It is agreed that one plotting sheet should be used for each 
reach, as shown in Fig. 9. The form used in Fig. 2 was presented for the 
purpose of illustrating the application of the proposed method. The method of 
plotting shown in Fig. 10 is an ingenious solution to the practical difficulty in- 
volved in plotting a graph to cover a range of 50 ft in elevation ata section. It 
has been found that the method of inclusion of overbank flow, proposed in 
Eqs. 23 and Table 11, is a commonly used expedient. Basically, this method 
corresponds to the use of the function," 


“Handbook of Hydraulics,” by H. W. King, McGraw-Hill Book Co., Inc., New York, N. Y., 1939, 
pp. 530-531. 
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Although some engineers may object to this procedure on the grounds that the 
value of A R} for a cross section taken as a whole is not equal to the sum of the 
(A R})-values of its component parts, it should be realized that the difference is 
not great for moderate amounts of overbank flow. If the value of n is deter- 
mined from observed flood profiles (using the sum of the (A R!)-values for 
channel and overbank), the computation of surface profiles should also use this 
n-value. 

Mr. Buehler is correct in his interpretation of the procedure to be followed 
in applying the proposed method to channels where it is necessary to use a 
different roughness coefficient for each reach. The computations in Table 
11 and Eq. 23 show how it is possible to use different roughness coefficients at 
the same cross section for different parts of the cross section. The writer agrees 
with Mr. Buehler that it is a cumbersome operation to determine the value of 
the roughness factor from flood observations when there is appreciable over- 
bank flow. 

Mr. Steinberg makes a welcome suggestion regarding the use of a weighted- 
area curve that can be applied directly to the discharges to obtain weighted 
velocities when overbank flow is involved. It has been observed, however, 
that the use of Eq. 5 involves less work and is simpler than the use of Eq. 9c or 
Eq. 10 to prepare a weighted-area curve. 

Mr. Steinberg’s discussion makes it obvious that it was not made sufficiently 
clear that the proposed method is a rigorous, although semigraphical, solution of 
the basic equation of varied flow given by Eq. 1. With this solution, it is not 
necessary to ignore the velocity-head terms or to make trial-and-error computa- 
tions. The method, although perfectly applicable to channels of uniform 
cross section and roughness, is especially suited to natural channels in which 
(1) the velocity varies from reach to reach, (2) consecutive lengths of reach are 
not equal, (3) each reach has a different roughness coefficient, (4) conditions of 
overbank flow prevail, (5) the same cross section has different roughness coeffi- 
cients for different water-surface elevations, and (6) the effect of eddy losses 
and bridge-pier losses can readily be included when necessary. Mr. Steinberg 
is correct in showing that, in most graphical and semigraphical methods which 
seek to eliminate trial-and-error computations, the velocity head corrections 
have not been taken into account. That is, these methods solve Eq. 1 by 
omitting consideration of the terms and es and furthermore by assuming 
that the difference in elevations Z, and Z, at the two ends of a reach are equal to 
the friction-head loss only. Although these simplifying assumptions are useful 
and not inaccurate when the velocity heads at each end of a reach are equal 
(in which case there is no varied flow), or when the difference between velocity 
heads is small compared to the friction-head loss in the reach (which is not the 
case when short lengths of reach are used under conditions of varied flow), they 
cannot be termed a rigorous attempt to solve the equation of varied flow. 

Mr. Steinberg presents a method of using a graphical solution for the 
equation of varied flow. Since the computation of M from Eq. 15a requires 
the determination of the properties of a midsection, as well as of the ends of a 
reach, the additional work involved is not worthwhile since no greater precision 
is obtained as Table 9 shows in comparison with Table 6. 
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DISCUSSION OF DAM MODIFICATIONS 
CHECKED BY HYDRAULIC MODELS 
PROCEEDINGS-SEPARATE 184 


G. H. Hickox,’ M. ASCE.—A praiseworthy description of an unusually 

7 Associate Director, Eng. Experiment Station, Univ. of Tennessee, Knoxville, Tenn. 
complete set of model tests has been presented by the authors. This paper 
presents an excellent example of the type of model testing that is necessary for 
the design of a new dam or the reconstruction of anolddam. The authors have 
arrived at completely satisfactory conclusions. Only a few comments, accord- 
ingly, on the general subject of model tests for dams shall be added. 

It has been assumed that considerable erosion would occur in the river bed 
below a dam, even in rock, if special provisions were not made to prevent it. 
There is considerable evidence to support this assumption. Wilson Dam on 
the Tennessee River, near Florence, Ala., is founded on limestone but it was 
not provided with an apron. Considerable erosion occurred at the toe when 
the dam was first placed in operation. Extension of the apron served merely 
to transfer the erosion to a point further downstream. It seems probable that 
properly conducted model tests might have abolished this problem. Messrs. 
Harrison and Kindsvater describe a dam that is twenty-five years old, which, 
although lacking an apron, has sustained considerable flow without any sig- 
nificant erosion of the river bed. The experiments were conducted on models 
having erodible beds in which considerable movement occurred. The results 
of the model tests were used for comparative purposes only, however, and no 
attempt was made to claim that erosion in the prototype would agree quanti- 
tatively with erosion in the model. The fundamental difficulty appears to be 
a lack of knowledge about the erodibility of various rocks under the flow condi- 
tions to which they are subjected. Without such knowledge, the designer 
must be cautious. 

The authors made tests on models at two different scales for various pur- 
poses. These objectives were to study the following: (1) Effect of the buttresses 
on erosive tendencies; (2) effect of sills at the end of the apron on erosion; (3) 
most suitable schedule of gate operation; (4) location of struts between but- 
tresses; (5) elevation of the powerhouse training wall; (6) elevation and loca- 
tion of construction cofferdams; (7) pressures on the face of the spillway; (8) 
water surface profiles; and (9) spillway discharge calibration. 

A large model of a section of the spillway was built at. a 1/40-scale to deter- 
mine the last three objectives. A smaller, 1/80-scale model of the entire dam, 
powerhouse, and a part of the river channel was built for the first six objectives. 
The paper indicates clearly that it is not sufficient to build a large-scale model 
for one or two bays of a spillway to study the erosion at the toe of a dam. 
There is usually a problem of erosion at the end of the spillway adjacent to 
the nonoverflow section of the dam or adjacent to the powerhouse. When- 
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ever the overflow section of the spillway is away from either bank, an eddy 
current is set up in the region of relatively quiet water between the end of the 
spillway and the adjacent river bank. At the base of the dam, or in the power- 
house tailrace, the eddy current moves toward the spillway section where the 
water surface is depressed. At the junction of the high velocity stream from 
the spillway and the lateral current from the eddy, local eddy currents of con- 
siderable intensity and erosive power are formed. These local eddies cause 
scouring at the end of the spillway or at the training wall (if a training wall 
exists at that point). Extension of the training wall in a downstream direction 
does not eliminate the difficulty except at a prohibitive cost. The solution to 
the problem usually requires a training wall with a proper combination of 
height, length, and location, and river bed paving. The proper combination 
can be determined only by a study of a complete model of the entire spillway 
and river and cannot be determined by tests made on a sectional model. 

The Bartlett’s Ferry Dam is unusual with respect to the relative locations 
of the powerhouse and the spillway. At this dam the powerhouse tailrace is 
directly below a part of the spillway. The training wall dividing the spillway 
from the tailrace is subjected to an increase in both water surface and water 
pressure on the spillway side, rather than reduction in elevation, as is usually 
the case when the dam is straight. 

The results reported as having occurred with the use of a sloping sill at the 
end of the apron are in agreement with the experience of the Tennessee Valley 
Authority (TVA). Aprons of most TVA dams are terminated with a sloping 
sill which reduces the erosive tendencies at the end of the apron. TVA experi- 
ence also agrees with that reported by Messrs. Harrison and Kindsvater that 
the slope of the upstream face is usually a compromise between the depth of 
erosion and the height of the waves produced in the river channel below the 
dam. The steeper sills produce the least erosion and the highest waves. The 
flat sills increase erosive tendencies but reduce the wave height. The slope 
adopted is generally chosen with regard to the importance of erosion at the toe 
of the sill and the effect of waves on navigation and erosion of the bank. 

It is encouraging to note that the authors considered conditions during 
construction. Model tests of high dams are frequently made without regard 
to conditions existing during diversion. A stilling basin that may be entirely 
satisfactory for the maximum flood discharge of a stream after the completion 
of the dam may be wholly inadequate when it is required to handle a much 
smaller flood under entirely different conditions. During the construction 
period, flow may be diverted through a narrowed river channel or, occasionally, 
through tunnels bypassing a cofferdam. In either case, the flow is concen- 
trated, and considerable damage may result. This is especially true when 
diversion occurs through spillway blocks which have not reached their final 
height. This circumstance causes concentration and direction of flow that 
are quite different from those of the completed structure, and the erosive 
effects may be quite unexpected. The concentration of flow over a fraction of 
a long apron may result in eddies that move debris and gravel onto the apron 
and, by continual movement of this material, cause serious scour. Such 
conditions can be predicted only by simulating construction conditions and 
expected floods with models. 
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It is gratifying to note that observations on the model have been made for 
rather small protype discharges. One of the reasons why satisfactory com- 
parisons of model and prototype operations are uncommon is that many models 
are tested only for the worst design condition, which is usually at or near maxi- 
mum discharge. Extreme discharges occur very rarely in the prototype. 
Comparisons of model and prototype are therefore available only when the 
experimenter has had sufficient foresight to conduct tests in the normal operat- 
ing range. The meager stock of knowledge about the reliability of models will 
increase if careful test data are obtained in the range of normal, or relatively 
frequent, prototype operating conditions. 

Another somewhat neglected phase of model testing is the effect of gate 
operations on erosive tendencies. As indicated by the authors, concentration 
of flow is usually undesirable and nearly uniform distribution is greatly pre- 
ferred. Tests on a model with an erodible bed are an extremely valuable 
method of determining the erosive effect caused by various methods of gate 
operations. Similar tests with a fixed-bed model are also valuable in determin- 
ing the effect of concentrated discharges on surface currents which may hinder 
or improve navigation. It would seem preferable to learn these combinations 
by model testing rather than by interfering with navigation. 

The importance of having model tests designed, carried out, and interpreted 
by a competent, well-trained staff should be emphasized. It is known that in 
any argument the conclusions reached are always a function of the initial 
premises. If the model experimenter assumes that only a single sectional 
model will be sufficient for the problem at hand, the answers obtained will 
be a function of the particular model. For example, such a model would 
be adequate for determining the effect of the height of the sill on erosion. If, 
however, it is assumed that there may be lateral flow problems, a complete 
model should be built to explore the situation. In such a case, the results dis- 
covered may be other than those predicted by use of the sectional model. 
Similar considerations apply to the use of fixed and movable beds. Model 
tests should therefore be made and interpreted by the most competent in- 
vestigators available. The requirements should include both field experience 
with operating structures and training in fundamental hydraulics and fluid 
mechanics. 


E. S. Harrison,* anp Cart E. KinpsvaTer,? Mempers, ASCE.—The 


® Asst. Chf. Engr., Georgia Power Co., Atlanta, Ga. 
* Prof. of Civ, Eng., Georgia Inst. of Technology, Atlanta, Ga. 


writers are indebted to Mr. Hickox for his comprehensive discussion of the 
model testing aspects of the paper. His general comments and admonitions 
are an invaluable complement to the paper. The writers also wish to acknowl- 
edge their use of the uplift measurements made at Bartlett’s Ferry Dam” 


10 “Uplift in Masonry Dams’’: Final Report of the Subcommittee on Uplift in Masonry Dams of the 
Committee on Masonry Dams of the Power Division, 1951, Transactions, ASCE, Vol. 117, 1952, p. 1243. 


(Georgia). 


| 
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DISCUSSION OF SETTLING RATES OF SUSPENSIONS 
IN SOLIDS CONTACT UNITS 
PROCEEDINGS-SEPARATE 186 


Norvat E. Anperson,* M. ASCE.—The generally accepted theory of 

6 Engr. of Treatment Plant Design, The Sanitary District of Chicago, Chicago, III. 
quiescent settling has been ably presented by Mr. Kalinske. The author also 
has presented mathematical functions for some of the variables affecting sedi- 
mentation (including upward flow) and has given some practical suggestions 
for good design (based on empirical results). However, he has not thoroughly 
presented general information on the settling characteristics of various sus- 
pensions, nor has he completely shown how to apply this data to the design and 
operation of water and waste treating and clarification basins. Data available 
for use in the equations are few. 

If, for example, a designer decided to use Ey. 10 to determine the settling 
area of an upflow basin, the only known factor would be Q, the inflow. The 
designer could assume C,, the solids concentration in the inflow, at 2,000 ppm, 
and C,, the solids concentration of the slurry, at 10,000 ppm. This makes the 
ratio C,/C, equal to 1/5. The subsidence rate V or V, must then be deter- 
mined. If samples of the suspension were available to make laboratory deter- 
minations of the rate of subsidence, there would still remain the question as to 
how the laboratory rate would be affected in the plant-size unit by stray 
currents, and how much of a factor should be allowed for bulking or variations 
in the behavior of the suspension. Hence, the designer may assume a surface 
rate of 1.5 gal per min per sq ft, or some other empirical figure based on ex- 
perience with similar suspensions. 

Data on the size, shape, density, and other characteristics of the mass of 
heterogeneous particles encountered in water and waste treatment are com- 
paratively meager. It is questionable whether there will ever be sufficient 
data of this kind to make a rigorous theoretical analysis of the design of a 
settling tank for a sewage plant based on the internal mechanics of suspensions. 
This is because the suspended particles in a single sample are so variable that 
they defy practical classification. The settling characteristics of these particles 
change greatly (sometimes within a few hours) in the same plant, and there 
are often considerable differences in the settling characteristics in different 
plants handling the same type of wastes. This does not mean that the theo- 
retical approach should be neglected. Knowledge of the mathematical re- 
lations, even though they cannot be evaluated numerically, provides a better 
understanding of empirical results, and sometimes lends new significance to 
old experimental data. 
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For example, W. Rudolphs and I. O. Lacy found that activated sludge 
settled more rapidly in a funnel than in a cylinder, and more rapidly in a 1-liter 
cylinder when set at an ang\e than when set vertica!.’ This could be explained 


7 “Settling and Compacting of Activated Sludge,” by W. Rudolphs and I. O. Lacy Sewage Works 
Journal, July, 1934, p. 647. 


by Eq. 8 and by the statement that 

“* * * Tn an expanding cone the faster settling particles will also remain 

towards the bottom but, beeause V’ and v both decrease upwards, this will 

tend to cause the ratio V’/v to change more gradually. Thus the suspension 
will tend to be of a somewhat more uniform concentration in a cone than in 

a cylinder.” 

Thus, if the concentration tends to remain more uniform, the heavier particles 
will continue to influence the settling rate of the lighter particles for a longer 
time, resulting in more rapid settling of the sludge mixture; or the observed 
differences in settling time could have been caused by the differences in the 
height of the vessels used in the test. 

Another example comes from the investigations of the behavior of the final 
settling tanks at the Southwest Sewage Treatment Works of the Sanitary 
District of Chicago, Ill, results of which were presented by the writer 
in 1944.% After it was discovered that density currents were the cause of floc 


®“Design of Final Settling Tanks for Activated Sludge,’ by Norval E. Anderson, Sewage Works 
Journal, January, 1945, p. 50. 


passing over the effluent weirs on the periphery of the 126-ft-diameter tanks, 
data from sludge-density determinations made in 1931 at the Des Plaines 
River plant in Illinois developed a new relation between the location of the 
effluent weirs and the sludge draw-off. 

With our present state of knowledge (1953), the approach to the design of 
settling basins through applied mechanies is of great value to the investigator 
who interprets results from operating units with a view toward improvements. 
The designer should be familiar with these fundamental relationships to apply 
intelligently the mass of empirical data which is available. A comprehensive 
discussion of the settling theory has been presented by Thomas R. Camp, 
M. ASCE. 


wails Satvontetion and Design of Settling Tanks,”’ by Thomas R. Camp, Transactions, ASCE, Vol. 111, 
p. 895. 


Mr. Kalinske is to be commended for his presentation of theoretical re- 
lationships, for emphasizing their importance, and for his suggestions for 
practical application. 

Ross E. McKinney,’ J. M. ASCE.—In an effort to present a compact 


- * Asst. Prof. of San. Eng., Dept. of Civ. and San. Eng., Massachusetts Inst. of Technology, Cambridge, 
Jlass. 


description of the action of the solids contact units, the author has tended to 
oversimplify the conditions encountered in such units. The descriptions of 
plain settling and upflow suspension were based on the results obtained with 
suspensions of sand, and it was concluded that the data could be applied to any 
type of suspension. Sand grains are dense, discrete particles which do not 
affect one another except by physical packing. The grains retain their identity 
as individual particles in concentrated suspensions (as they also do in dilute 
suspensions) and are unaffected by turbulence or any shearing action which 
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might be created by the velocity of water passing through the suspension. 
Unfortunately, sand does not represent the type of solids normally removed by 
solids contact units. 

In the clarification of water, the solids contact units are used to remove 
finely divided suspended matter and colloidal matter with the aid of chemical 
coagulation. The floc that is formed is not a tough, discrete particle, such as 
a grain of sand, but is a fragile particle easily affected by other particles and by 
liquid motion. Chemical floc or biological floe produced during the clarification 
of sewage is of the same nature. In solids contact units, the upflow rate is 
not only a factor of the settling velocity of the particle but also a factor of 
the structural strength of the floc. 

At low rates of flow, the quantity of floc formed is small and the upflow 
velocity is low. Thus the denser flocs will settle readily, leaving a small 
concentration of light floc in suspension. As the flow increases, the quantity 
of fioc and the upflow velocity also increase. The heavier floc is retained in 
suspension, and—since more floc is present—the number of contacts with other 
floc particles increases. If the upflow velocity is less than the shearing velocity 
of the floc, larger particles will form. The larger particles (increasing in mass 
and size) will oppose the upflow velocity and tend to become more compact. 
The suspended floc will form a “blanket” if the floc concentration is sufficiently 
high. The height at which this blanket will remain suspended in the unit 
will depend on the velocity of flow through the blanket without damaging the 
structural integrity of the blanket. The fragility of the floc is such that it is 
not possible in normal operation (without some destruction of the floc that 
forms the blanket and a resultant decrease in efficiency of operation) to increase 
the upflow velocity to the point where the blanket would be lifted to the surface. 

Eqs. 9 and 10 are valid only if the units are properly chosen. The total 
weight of the settling solids, expressed in Eq. 9, has two variable factors, A 
and V. The settling area is usually taken as the surface area of the settling 
section of the solids contact unit. This is correct only if the solids are dis- 
tributed uniformly over this area which will be a factor of the inlet design. 
The settling velocity of the floc is the subsidence rate of the sludge blanket. 
* This velocity is affected by the floc structure, which in turn is affected by the 
temperature of water, the concentration of solids to be removed, the concentra- 
tion of chemicals added, the mineral content of the water, and the pH of the 
water. As a result of the nature of the floc, plain settling curves will not give 
the true settling velocity for Eq. 9. 

Transformation of Eqs. 9 and 10 into actual designs involves many practical 
problems. As stated by Mr. Kalinske, the type of solids contact unit shown in 
Fig. 5 is dependent on the flow through the unit for maintaining the solids level. 
The unit shown in Fig. 6 overcomes this by recirculating the slurry. This 
type of unit affords the additional advantage of introducing the solids into the 
settling unit with a definite downward velocity and thus utilizing the density 
currents to aid in settling. However, the unit shown in Fig. 6 has the dis- 
advantage of not distributing the solids uniformly over the entire settling area. 
This disadvantage is minimized by increasing the depth at which the solids 
enter the settling section. Thus, at best, the solids contact units will be (like 
most units in water and sewage disposal systems) a compromise. 
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The author has attacked rectangular secimentation basins as being of poor 
design. Mr. Kalinske suggests increasing the efficiency by increasing the depth 
of the sedimentation basin. This is contrary to the ideas of Mr. Camp, which 
have been incorporated into the water treatment plant at Cambridge, Mass. 
At Cambridge a conventional sedimentation basin was used to increase the 
sedimentation efficiency by adding two trays at 4 and j of the original depth. 
Thus the sedimentation efficiency was increased by reducing the depth rather 
than by increasing it. 

On the basis-of empirical design (and most water and sewage clarification 
units are so designed) there is a definite need for further research on the 
mechanisms involved. Sound theoretical concepts combined with practical 
design and operation will eventually yield the most efficient solids contact 
units as well as other types of solids separation units. 

A. A. Kauinske,"” M. ASCE.—The writer has attempted to indicate some 

10 Director of Research and Development, Infileo Inc., Tucson, Ariz. 
of the basic physical problems entering into the design of solids-contact basins. 
Messrs. Anderson and McKinney both state that the general relationships 
presented by the writer require much data regarding the settling characteristics 
of suspensions before an engineer can use such relationships in his design work. 
These data were not presented because they are nonexistent. Sych data are 
dependent (as Mr. McKinney mentions) on so many other factors that mere 
tabulation for one set of conditions would serve no useful purpose. It was 
shown that such data should be obtained directly by relatively simple labora- 
tory tests for the specific type of flocculated suspension in the basin. Of 
course, if it is impractical to make such measurements directly, the designer 
must examine data obtained from similar types of suspensions. 

Mr. McKinney criticizes the use of the theory of the settlement of suspensions 
of discrete particles to analyze the settlement of suspensions of flocculated 
particles. As both Messrs. Anderson and McKinney observe, the behavior 
of flocculated particles is complex; any attempt to formulate rational mathe- 
matical expressions for such behavior would unquestionably be useless for 
design purposes. When one encounters such complex phenomena, many sim- 
plifications must be made in any theoretical analysis. Such analyses’are none- — 
theless useful as a starting point for systematizing empirical data. 

It was carefully emphasized that accurate design information could be 
obtained best from laboratory settling tests of suspensions under conditions 
approximating those (for example) in an upflow treating unit. Hundreds of 
such subsidence tests have been performed on a wide variety of flocculated 
suspensions of different concentrations. Such data have been extremely useful 
in the design and the capacity rating of solids-contact units. Subsequent 
study of the operation of such units in the field has confirmed the basic sound- 
ness and practicability of using such data as design criteria. 

Mr. McKinney states that “at low rates of flow, the quantity of floc formed 
is small,” and adds that therefore the quantity of floc present in the suspension 
will be low. In a properly designed and operated solids-contact treating basin, 
the total quantity of floc present in the blanket or slurry pool should be prac- 
tically independent of the flow rate. It is desirable (and possible) to keep the 
floc suspension at approximately a uniform concentration, irrespective of the 
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flow rate or the rate of floc production. This is accomplished by wasting solids 
from the suspension in proportion to their rate of production ; in many instances, 
this is in proportion to the flow. 

The operation of the solids-contact unit shown in Fig. 6 has been noted 
by Mr. McKinney, who comments on the fact that the downward-circulating 
slurry is aided in settling by a density-current effect. There is no settling of 
the solids in the circulating zone of this unit; this occurs only in the concen- 
trators, which occupy a part of the clarification area. There is a density- 
current effect, however, which helps in keeping the moving slurry from intermix- 
ing with the lighter, clear liquid above. A solids-contact unit, so designed, does 
distribute the solids uniformly over the entire clarification area. The circu- 
lation of the slurry at a rate equal to several times the treating flow rate is 
effective in distributing the solids equally over the entire clarification area of 
the basin. 

Mr. McKinney states that the writer “has attacked rectangular sedimen- 
tation basins as being of poor design.” This is not completely accurate. The 
writer did indicate that rectangular basins without proper take-offs may require 
larger and deeper basins to compensate for poor outlet design. Regarding the 
use of ‘‘trays’”’ in rectangular basins to increase their capacity, R. Eliassen, 
M. ASCE, in his discussion of Mr. Camp’s paper, pointed out some of the 
practical operating and design problems in connection with such design. It 
should be realized that this paper was concerned only with the design of solids- 
contact basins, which by their nature operate differently from conventional, 
straight, flow-through basins. 
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DISCUSSION OF SNOW HYDROLOGY FOR 
MULTIPLE- PURPOSE RESERVOIRS 
PROCEEDINGS-SEPARATE 189 


Ray K. Linstey,'® 4.M. ASCE.—The accumulation of snow in the moun- 

18 Associate Prof. of Hydr. Eng., Stanford Univ., Stanford, Calif. 
tains of the western United States is economically important. Without this 
snow the development of this area would have taken a far different and prob- 
ably less favorable course. To the engineer, however, this snow constitutes a 
troublesome design factor. Many approximations and assumptions are neces- 
sary in any computation of the runoff rates to be expected from melting snow. 
Mr. Riesbol has performed a service in presenting a fine summary of snow- 
hydrology practices. — 

The hydrothermogram procedure outlined by the author (under the head- 
ing, ‘‘Application to Spillway Design: Inflow Design Flood Based on Seasonal 
Snow Melt”’) is an interesting commentary on the present (1953) status of snow 
hydrology. The method is an application of the ‘“degree-day” approach in 
which temperature above some base value is taken as an index of the melting 
rate. It is similar in many respects to the procedure suggested by the writer'® 


16 “A Simple Procedure for the Day-to-Day Forecasting of Runoff from Snowmelt,"’ by R. K. Linsley, 
Transactions, Am. Geophysical Union, Pt. 3, 1943, pp. 62-67. 


and basically the same as that proposed by G. D. Clyde,!’?7 M. ASCE. 


17 “Snow-Melting Characteristics,"’ by G. D. Clyde, Bulletin No. 231, Utah Agri. Experiment Station, 
Logan, Utah, 1931. 


A great deal of fundamental research in the field of snow physics has been 
conducted. Many features of the subject are imperfectly understood but a 
skeleton theory exists which should permit one to go beyond the degree-day 
approach to include in the melting index such parameters as solar radiation, 
humidity, and wind. 

Mr. Barnes did not use these parameters (see under the heading, ‘‘Applica- 
tion to Spillway Design: Inflow Design Flood Based on Seasonal Snow Melt’’) 
because the necessary data for their use on an engineering scale were not avail- 
able to him. A region of heavy snow accumulation is a region of scanty meteor- 
ological and hydrologic data. The variability of meteorological and hydro- 
logie elements in rugged terrain and the economic importance of the region 
suggest that a station density greater than the national average might be in 
order. Continued research will fill out the skeleton of knowledge regarding 
snow physics. With this knowledge the basic data required to utilize it should 
also be available. 

The amount of snow melted from a given basin in any period of time depends 
on the area of the snow pack exposed to melting conditions. This area is 
bounded by the lower edge of the snow cover and by the upper limit of the 
melting temperatures. Neither of these boundaries is strictly an elevation con- 
tour. Meteorological data may suffice to- define the level to which melting 
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temperatures extend, but the lower limit of the snow cover can be determined 
only by field location. Aerial photography offers one of the best means of 
obtaining these data, although terrestrial photography and reconnaisance map- 
ping methods contribute to local determinations. It seems axiomatic that a 
rational computation of snow melt cannot be performed without some knowl- 
edge of the controlling area of the melting zone, just as the runoff from rainfall 
on a basin cannot be estimated without first knowing the drainage area. 

It is necessary in some cases to vary both the base temperature for the 
hydrothermogram and the degree-day factor, or melting rate of runoff, as the 
season progresses. The variation in the base temperature results from the 
fact that the melting zone moves progressively higher above the temperature 
station as the melting proceeds. Similarly, the variation in the melting rate 
of runoff must result from the change in area of the melting zone. Other 
factors such as the seasonal variation in solar radiation are undoubtedly impor- 
tant, but there is clearly a need for data defining the critical melting zone. 

In Fig. 1 is shown an analysis which has been applied to the computation 
of the snow-melt hydrograph of a small area. When this analysis is applied 
to the hydrographs of snow-melt runoff from large areas, it is often found that 
the recessions of the daily increments of snow melt intersect the falling limb of 
the total hydrograph. A more rational approach seems to visualize the total 
hydrograph as the sum of a hydrograph of direct runoff and a ground-water 
(or subsurface flow) hydrograph (Fig. 6(a)). This approach makes apparent 
the many problems which hamper the application of experience gained in the 
study of small basins to the engineering problems of large basins. It is quite 
probable that the development of techniques must wait for available data on 
large basins. 

Under the heading, ‘‘Capacities of Multiple-Purpose Reservoirs: Runoff 
Forecasting Curves,’”’ Mr. Riesbo! mentions only the use of snow surveys, 
whereas successful forecasts of seasonal runoff have been made using precipi- 
tation data recorded at regular stations of the Weather Bureau.'® Existing 


18**Recent Developments in Water Supply Forecasting from Precipitation,’’ by M. A. Kohler and 
R. K. Linsley, Transactions, Am. Geophysical Union, 1949, pp. 427-436. 


records of rainfall make possible forecasts of runoff for areas where snow 
surveys are not available. These records offer a better source of data for 
statistical analysis than do the snow-survey records. Since the data on pre- 
cipitation are regularly available, the cost of the forecasts is less than the cost 
of forecasts made from snow surveys. It appears that the accuracy of the 
forecasts based on precipitation data is as good as those based on snow surveys, 
although the superiority of either method depends on the local area and the 
data available. Precipitation data for forecasts of this type certainly deserve 
consideration for use in snow hydrology. 

R. W. Gerpvet," M. ASCE.—An improved approach to the development 


1 Chf., Field Operations and Analysis Branch, Snow, Ice and Permafrost Research Establishment, 
U.S. Army, Wilmette, II. 


of operational procedures for multiple-purpose reservoirs in areas where 
precipitation in the form of snow is the dominant hydrologic element has been 
presented by Mr. Riesbol. The author has emphasized the fact that the 
hydrologic engineer is compelled to use empirical relationships to develop 
suitable design and operation criteria despite the considerable amount of 
research that has been conducted on snow-melt and runoff problems. The 
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gaps in the scientific literature concerning snow may be sufficiently large to 
justify the use of empirical relationships, but the use of hypotheses not com- 
patible with experimental evidence or the known physical properties of snow 
is not justified. 

The literature on snow hydrology, mechanics, and physics contains sufficient 
references to the liquid-water holding capacity of snow and to the metamorphic 
processes which occur in snow to refute any assumption that a snow pack must 
achieve a density of 40% before discharge of rain or snow-melt water will take 
place. Likewise, there is no basis for assuming that a snow pack will not 
continue to increase to a density of more than 40% without storage of liquid, 
or that it will not retain that increase in the form of the solid phase. 

In 1948 the writer examined the existing pertinent literature in a search for 
information about the physical properties of a ripe, discharging snow pack, for 
a report to the International Commission on Snow and Glaciers.” It was found 

2» ‘Physical Changes in Snow-Cover Leading to Runoff, Especially to Floods,"’ by R. W. Gerdel, 
Union géodesique et géophysique internationale Association d'Hydrologie Scientifique, Assemblée gén- 
érale d'Oslo, Vol. II, Oslo, N. rway, 1948, pp. 42-54. 
that densities between 25% and 56.5% for a ripe, melting, actively discharging 
snow pack were reported by several investigators. A snow pack usually consists 
of a number of readily identifiable layers which are the products of separate 
storms. The variation in density between these layers may be large at the time 
of deposition and will be maintained even in the ripe pack. R. A. Work found 
that homogeneous density was never obtained in a ripe pack and that the 
density of different layers in a late-April, ripe snow pack varied between 45% 
and 58%." In the Central Sierra Snow Laboratory observers found variations 


"Snow Layer Density Changes,” by R. A. Work, Transactions, Am. Geophysical Union, Vol. 29, 
194, pp. 525-546. 


of more than 10% between the densities of ripe, discharging, spring snow packs 
at different sample points in the 4-sq-mile basin. 

The results of investigations of the discharge of water from a melting snow 
pack made by the Weather Bureau at the Soda Springs Cooperative Research 
Project, in California, showed that a pack of 52.8% density, rapidly melting 
under an air temperature of 65°F, had a total liquid and solid water-equivalent 
content of 14.8 in. at 3:30 p.m.” This content consisted of 13.4 in. of water- 


22 “Snow Studies at the Cooperative Snow Research Project, Soda Springs, California,’’ Annual Report, 
U.S.W.B., 1943-1944, pp. 11-14. 


equivalent snow crystals and 1.4 in. of liquid water. At 8:30 the following morn- 
ing, when temperatures had been low enough during the night to prevent or re- 
tard melting, the pack density was still 53%, whereas the total water equivalent 
was only 13.8 in., consisting of 13.4 in. of water-equivalent snow crystals and 0.4 
in. of liquid water. The studies showed that no liquid water was frozen in the 
snow during the night and that the 1-in. loss in water equivalent of the pack 
during the night could be accounted for by the drainage of the 1 in. of liquid 
water from the pack during the same period. The maintenance of almost 
constant density was due to shrinkage and reduction in voids. 

The 1944-1945 report from this project presented the results of investiga- 
tions of a heavy rain on new and old snow. A new snow, five inches deep, and 


2 “Snow Studies at the Cooperative Snow Research Project, Soda Springs, California,’’ Annual Report, 
U.S.W.B., 1944-1945, pp. 10-13. 


of 10% density had been deposited on top of an old pack of 38.3% density. In 
places, the new snow was deposited on clean black top and concrete roadways. 
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During one day, 1.44 in. of rain fell upon the new snow. At the height of the 
storm the new snow, overlying the old snow, contained only 13% of liquid water. 
On the paved roadway, where water was impounded in depressions, the slush 
that developed had a free-water content of 30%. Under no circumstances 
could such a high free-water content be attributed to the water-holding capacity 
of thesnow. It was definitely the result of snow floating on, or being submerged 
in, water impounded in a depression with an impervious bottom. It is of 
interest to note that when efforts were made to measure the density of the 
slush by the usual tube sampling and weighing procedure, the water drained 
out of the sample so rapidly that the maximum measured density was only 
16.5°,—much less than the 40% density assumed to be required before melt or 
rain water could be discharged from the pack. + 

During the period between 1943 and 1950, when the writer was actively 
engaged in snow-hydrology research in the western United States for the 
Weather Bureau and the Corps of Engineers, a number of experiments were 
performed to determine the effect of large amounts of water on the shrinkage 
and increase in density of monolithic cylindrical samples of snow. A cylindrical 
container equipped with a cutter was used to obtain snow core samples (8 in. 
in diameter and 24 in. long) from sensibly homogeneous horizons in the pack. 
A screen bottom was applied to the container, the sample was weighed, and 
was then saturated with water at 32°F. The core was returned to the pack to 
establish natural capillary tension through the screen bottom. At the end of 
approximately 4 hr, the core was removed from the container, carefully 
calipered and weighed, and the free-water content was measured by the 
calorimetric method. These studies showed that an 18%-density dry snow at a 
temperature of 32°F underwent a reduction in volume and thus increased its 
density to 26% with an unfrozen water content of 1.5%. During the shrinkage 
process more than 2 in. of the previously applied water drained out of the 
sample, but there was no loss of solids. A 24%-density snow was increased to 
a density of 31% by this process. Repeated addition of liquid water, equal 
to three times the water equivalent of the original solid snow crystals in the 
core sample, changed the density from 22% to 34% by reducing the volume 
while allowing more than 15 in. of liquid water to pass through the initial 
24-in.-long core of snow. The water equivalent of the sample at the beginning 
and end of the experiment was 5.28 in. and 5.48 in., respectively. The positive 
increase of 0.2 in. in water equivalent probably represents the amount of liquid 
retained in the core. 

The results of these studies have not been published heretofore because the 
investigations were limited to high mountain snow covers where steeply 
sloping land provided rapid drainage and because there was lacking the time 
and equipment to collect sufficient data with which to demonstrate the extent 
of application of the results. They are presented in this discussion in hope 
that the 40%-hypothesis may be subject to a revision based on experimental 
evidence rather than on empirical judgment. 

In another paper™ the writer has discussed some experiments performed on 


“The Transmission of Water Through Snow,” by R. W. Gerdel, Minnesota Eepprmationst Hydraulics 
Gusdain Univ. of Minnesota, Minneapolis, Minn., August 3i -September 4, 19. 


the water transmission and retention capacity of snow as eee with an 
electronic capacitance meter developed for this type of investigation. It was 
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found that for snow of 35% to 46% density, water applied to the surface of 
the snow (in quantities of 0.5 in. to 2.0 in.) was transmitted througi: 0.9 in. to 
24 in. of snow per min. This velocity appears to indicate a limited temporary 
detention capacity. The permanent retention capacity, when free drainage was 
available, was approximately 0.1 in. of water per ft of 50%-density snow. It 
was shown that the movement of water through the snow pack follows definite 
channels which develop perpendicularly and horizontally within the pack. 
These channels have been recognized by glaciologists who term them “firn 
pipes” or “glands.” The channels develop into drains filled with coarse-grained 
crystals in the melting seasonal layer of snow and are similar to the rock-filled 
drains used in agricultural and other soil-engineering practices. When the 
ground is exposed as a melting snow pack recedes in the mountains, the flow 
channels which have developed at the contact between the snow pack and the 
ground are evident on some soils, as a dendritic pattern of erosion marks. 
When these channels develop over rock or less erodible soils, they appear as 


arched piping in the bottom snow layer. 
There is little experimental evidence to indicate that a deep and dry snow 


cover will retain appreciable amounts of rainfall, as stated by the author. 
The presence of free water induces a rapid metamorphic process which reduces 
the surface area of the individual snow crystals to the minimum area compatible 
with the mass and stable form of the crystal. That is, the feather-like, or 
dendritic, snow flakes are converted to small crystals of hexagonal form which 
resemble a mass of table salt or rock salt. The surface detention capacity of 
the snow flake is rapidly reduced by the metamorphic process, and the water 
which may have been retained on the larger surface area is discharged. This 
metamorphic process is not uniform throughout the pack. The process 
appears to proceed more rapidly at one place than at another, a phenomenon 
which accounts for the development of the flow zones within the pack and at 
the ground surface. 

The size of the metamorphosed grains in a spring snow pack is the same as 
the size of a mixture of fine and coarse sands, about 0.02 mm to 2.0 mm in 
diameter. A reference to the literature on soil moisture shows that a sand 
with such a size distribution and free drainage will hold from 1% to 3% of liquid 
water, the same quantity as has been measured and reported for a spring snow 
pack. 
The temporary detention of liquid water in snow is a function of the time 
required for the development of internal flow channels; hydrologically, these 
channels are effective temporary supplemental tributaries to the normal 
drainage pattern of the basin. There is no evidence that snow must achieve 
a 40% density before this supplemental tributary system becomes operative 
nor that, having achieved this density, the supplemental tributary drainage 
can be considered as having a discharge capacity which is a function of the 
density in excess of 40%. 

In presenting this discussion the writer does not wish to detract from the 
value of Mr. Riesbol’s paper. This discussion is submitted so that hydrologists 
may give more recognition to the known physical and thermodynamic properties 
of snow and the -associated hydraulie phenomena in developing snow-melt 
runoff formulas and design and operational criteria for multiple-purpose 


reservoirs. 
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As a source of information, attention is called to a bibliography of six: 
thousand items on snow, ice, and frozen ground available in the Engineering 
Societies’ Library, New York, N. Y.2° 


2? ‘Annotated Bibliography on Snow, Ice and Permafrost,” SIPRE Report 12, Vol. UI, Snow, Ice and 
Permafrost Research Establishment, Corps of Engrs., U.S. Army. Prepared by the Library of Congress, 
Washington, D. C., January, 1953. 


Harotp D. Harterson,?* A.M. ASCE.—One of the most important as- 


2% Chf., Hydrology Branch, Project Planning Div., Bureau of Reclamation, U.S. Dept. of the Interior, 
Region 1, Boise, Idaho. 


pects of the multiple-purpose use of major reservoirs in the western and north- 
ern sections of the United States has been investigated by Mr. Riesbol. This 
subject was previously examined in a symposium—one of the papers of which 
was written by Albert L. Cochran,?7 M. ASCE. 


2 “Their Use for Flood Control,” by Albert L. Cochran, from ‘‘Multiple-Purpose Reservoirs; A 
Symposium,”’ Transactions, ASCE, Vol. 115, 1950, p. 818. 


The ability to forecast the volume and shape of the flood-runoff hydro- 
graph is a valuable asset in operating multiple-purpose reservoirs and in 
realizing the purposes which these reservoirs have been designed to meet. The 
author rightly states (under the heading, ‘‘Conclusions’’) that 


“Great refinement in procedure, with corresponding economy in both design 
and operation, will result from * * * the extension of networks for the 
collection of data so that they include pertinent elements on a watershed 
basis.” 


Continual study of forecasting procedures and operating criteria—to utilize 
available information on snow depths, areal coverage, and. melting processes— 
is being made (1953) by all reservoir-operating agencies in areas where snow- 
melt runoff provides a large part of the storable flow. Snow surveys are con- 
stantly being conducted to cover ever-widening areas. New methods of 
measuring depths and water content of the snow pack are being devised. Ex- 
periments are being continued to determine the relationship of snow melt to 
some of the more easily measured meteorological factors such as temperature, 
air movement, and relative humidity. 

It is worthwhile to note that the development of forecasting procedures, 
and the application of resulting forecasts to the operation of multiple-purpose 
reservoirs, have provided methods whereby relatively full use of certain 
storage facilities can be made for functions which are ordinarily in opposition 
to each other. For example, in certain areas it is possible to use the full 
capacity of certain reservoirs both for flood control and irrigation, or for 
other conservation purposes, despite the fact that flood control requires reserva- 
tion of reservoir space whereas irrigation demands the assured accumulation of 
storage for subsequent use. 

Mr. Riesbol’s use of the term ‘“‘operation ‘parameters’ ” might readily have 
been clarified by a definition of exactly what these parameters denote. The 
mathematical definition of a parameter is any constant or element whose values 
characterize one or more of the variables entering into a system of expressions 
or functions. Thus, the operating curves (which the author referred to under 
the heading, ‘Capacities of Multiple-Purpose Reservoirs’) show the relation- 
ship between the storage space required to control the forecasted flood within 
specified limits of outflow and the volume of runoff which is anticipated be- 
tween any given dates and the normal end of the flood season. The curves 
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shown in Fig. 2, empirically derived from floods of record, are curves that 
envelop the storage requirements for various volumes of total forecast runoff 
from any given date to July 31. An additional safety factor has been added 
to Fig. 2 to assure adequate reservation of space at the beginning of the flood 
period. This factor was added to account for the possible forecasting errors 
shown in the computations in Fig. 3. 

To be more complete, Table 1 might have included a split period beginning 
on April 9 and continuing through April 16 so that the negative values could be 
shown to balance the positive values of required storage space. By means of 
this inclusion, the summations shown in Col. 5, Table 1, would be correctly 
referenced as to date. 

The author’s presentation of three methods for applying snow hydrology 
to spillway design includes the best of the available (1953) methods for deriving 
spillway-design floods of snow melt or combination rain and snow-melt types. 
It is gratifying to note that Mr. Riesbol suggests either the use of the unit- 
graph method or the augmentation of an actual snow-melt event to distribute 
the design snow-melt volume. Both procedures are practical and are superior 
to the widely used scheme of proportioning runoff directly to a design thermo- 
gram, ignoring the runoff characteristics of the basin. 

It is felt that the detailed approach described in the second method (under 
the heading, ‘‘Application to Spillway Design: Flood Produced by Rain on 
Snow’’) is an admirable attempt to account for causative snow-melt factors 
other than air temperature. However, the number of assumptions necessary 
and the type of data required in the method prevent its practical use in most 
localities. Almost none of the numerous watersheds in the northwestern parts 
of the United States (for which spillway floods have been derived or will be 
derived by the USBR) has a temperature record within the watershed. The 
determination of values for K, U, and p in Eq. 3 and the determination of the 
required infiltration rates for such watersheds are highly speculative. 

It is advantageous to use the hydrothermogram method for the derivation 
of spillway floods for large watersheds where seasonal snow-melt floods of 
several months’ duration are prevalent. An examination of the method or 
methods for converting an actual and a design thermogram into corresponding 
runoff (through the hydrothermogram approach) would have been helpful to 
the potential user of this approach. 

In determining the relationship between temperature and snow-melt 
runoff (Table 2), a correction should be made for the area contributing to 
snow melt during a given recorded event. In some cases, snow-course data 
can be used to correct continuously the contributing area during the given 
snow-melt flood. The author states (under the heading, ‘Application to 
Spillway Design: Flood from Snow Melt Alone’’) 


“If the stream-gaging station is a considerable distance from the snow-melt 
area, it is nec essary to determine time lag between snow-melting tempera- 
tures and runoff.” 


The location of the stream gage is only one of the factors which can cause a 
lag between temperature and runoff. The location of the temperature station 
with respect to the snow-melt area is another factor. A third factor is the 
time which is required for the heat in the air to convert the snow into water. 
The first two factors cause variance in the apparent time lag observed in 
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different snow-melt floods because of differences in the direction and speed of 
the respective causative warm-air mass. The third factor varies, depending 
on the depth and quality of the snow. 

Should not the vertical axis of Fig. 4 be labeled, “Maximum fifteen-day 
melting rate of runoff, in inches per degree day” rather than “Maximum * * * 
of runoff, in inches?” 

In the application of the detailed and somewhat theoretical method for 
deriving a design flood produced by rain on snow (Table 4), it appears that 
consideration should be given to the melting of the snow by the heat of the 
rain water. This amount of melting can be computed from 


_ P (Te — 32) 
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in which M is the snow melt in inches, P denotes the precipitation in inches, 
and 7, is the wet-bulb temperature (assumed to be the temperature of the 
rain). Although generally regarded as insignificant, the heat in 7.2 in. of 
precipitation of the maximum 6-hr period at an assumed 50° F could cause 
0.90 in. of additional snow melt in the example cited in Table 4—thus increasing 
the contribution to the mean-water excess of the basin from 4.79 in. to 5.60 in. 
for that single 6-hr period. 

The relationship of snow melt to runoff has not been given as much con- 
sideration in the design and operation of multiple-purpose projects as its im- 
portance warrants. The qualitative use of available knowledge on this sub- 
ject, applied through the judgment of the engineer, should be bolstered by 
continuous improvement of procedures to apply this increasing knowledge 
quantitatively. 

Davip M. Rockwoop,”* A. M. ASCE.—The problems of snow hydrology 


_ % Technical Director, Snow Investigations, Corps of Engrs., U.S. Dept. of the Army, North Pacific 
Div., Portland, Ore. 


as related to the design and operation of water-control projects have been 
described by the author. Mr. Riesbol has also summarized methods for deter- 
mining rates of snow melt, reservoir capacities, seasonal forecasts, design floods, 
and storage of liquid water in the snow pack. 

The purpose of reservoir storage in the regulation of streams is twofold. 
This storage can be used for supplementing deficiencies in the natural flow for 
specified purposes, or it can be used for reducing overbank flows to specified 
channel capacities. In rivers where the eccurrences of low and high flows 
invariably follow a seasonal pattern, the reservoir-storage space can be used to 
meet both purposes compatibly, provided proper operational procedures are 
established. In rivers where the stream flow is primarily the result of snow 
melt, and when there is normally a delay between the accumulation and melt 
periods, optimum use of storage space can be made because the volume of 
runoff during the melt period can be forecast two or three months in advance. 
The Columbia River in the Pacific Northwest is a notable example of a stream 
with these characteristics and is particularly susceptible to the development of 
multiple-purpose reservoir storage. 

Efficiency of design and operation of multiple-purpose reservoir projects on 
snow-fed streams is dependent on the accuracy with which the stream flows can 
be forecast. In the matter of project design, the application of snow hydrology 
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involves determination of reservoir-storage requirements and the derivation 
of design floods to meet the criteria for spillways and appurtenances for the 
structure. For the operation of the project, knowledge of snow hydrology is a 
requirement for forecasting rates of flow for day-to-day operation, and also for 
volumetric-runoff forecasts in order to establish long-term operating criteria 
on a seasonal basis. 

Engineers of the Corps of Engineers have long recognized the importance 
of snow hydrology for efficient design and operation of their projects, par- 
ticularly those projects in the western United States. In 1945 the Corps of 
Engineers entered into a cooperative snow-investigation program with the 
Weather Bureau—with limited cooperation from other governmental agencies— 
for the purpose of determining criteria to be used in snow-hydrology problems. 
That program included (1) operation of three basin snow laboratories in the 
mountains of the western United States for periods of from 4 yr to 7 yr ,(2) 
processing and publishing the hydrometeorological data from the laboratories, 
(3) maintaining an analytical unit for performing basic research in snow 
hydrology, and (4) disseminating information to participating offices. The 
cooperative phase of the progratn ended in 1952; however, the Snow Investi- 
gations Unit was continued by the Corps of Engineers with offices located 
(1953) in the North Pacific Division office. The activities of the unit were 
directed, almost entirely, toward specific analytical problems which were divided 
into two categories. The first was the analysis of physical processes affecting 
snow deposition and snow melt. The second was the practical application of 
procedures to basin snow-hydrology problems, which dealt with forecasting 
rates of stream flow and estimating seasonal volume of runoff. 

The use of degree-day factors provides a quick method of approximating 
daily increments of snow melt from readily available meteorological data 
because air temperature is an index of heat transfer by several processes. The 
complexity and variability of these processes require that temperature-index 
methods be used with considerable judgment. This judgment can be obtained 
only by study and evaluation of the components of heat transfer for various 
meteorological and basin conditions. For example, on clear days the processes 
of heat transfer in a forest-covered area are quite different from those processes 
for open areas whereas the total daily energy exchange may be approximately 
equal. An air-temperature index in the forest area for this condition is much 
more reliable than for an open area because the air temperature is a less reliable 
index for solar radiation than for convective heat transfer and long-wave radia- 
tion from the forest canopy. Differences in meteorological factors, such as 
vertical and horizontal distribution of moisture and temperature of the air 
(together with the associated cloud layers, air drainage, and stability), are 
related to the general atmospheric circulation, which in turn produces changes 
in these properties with respect to time. All these properties affect the rate of 
heat exchange in a manner which cannot be defined by air temperature at only 
a few surface index stations. The problem of accuracy becomes a matter of 
economic necessity. The project design or operation must warrant the expen- 
diture of time and effoit to obtain the required basic data and to analyze the 
meteorological and basin effects for each day’s melt. When the requirements 
are such that empirical day-degree relationships are adequate, application and 
limitations of the method should be understood. 
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The examination by Mr. Riesbol of spillway-design floods resulting from 
snow melt alone is based on an analysis of historical stream-flow and tempera- 
ture records for a particular stream. Empirical melt rates were determined, 
anc interpolation to design conditions’ was made on the basis of a maximum 
observed runoff and a maximum enveloping melt rate corresponding to that 
runoff. These values do not necessarily represent the maximum flood poten- 
tial of the basin, when one considers such important factors as the maximum 
possible snow accumulation, the minimum basin losses, the most optimum flood- 
producing temperature sequence, the meteorological events leading to the 
greatest rate of heat transfer at the time of the flood peak, and the maximum 
snow cover that can reasonably be expected at the time of the peak. The 
simultaneous occurrence of these important flood-producing factors would 
probably lead to a much greater flood flow than that derived from the analysis 
of a relatively few years of stream-flow records—during which time it is quite 
unlikely that maximum flood-producing conditions have occurred. 

Floods produced by rain on snow are more complex and less certain of 
analysis’ than those caused by snow melt alone. Eq. 3 is used to evaluate 
heat transfer by convection and condensation but does”"not account for 
heat transfer ‘by solar or long-wave radiation. The application of Eq. 3 
to basin analysis depends on an evaluation of the basin constant, K. In 
the example shown for the Sly Park Dam site, it was not explained how K 
was determined for this particular drainage basin. The use of Eq. 3 for 
this example is valid, inasmuch as short-wave radiation is negligible under 
these conditions, and net long-wave exchange is primarily a function of 
air temperature. However, it should be realized that, for the evaluation of 
total heat exchange, radiation ordinarily causes a large part of the total melt. 

The retention of liquid water in the snow pack has been the subject of much 
speculation. Qualitatively, the pack has been considered a “sponge” 
capable of retaining large amounts of liquid water entering as a result of rain- 
fall or surface melt, thereby reducing the rates of flow that would have 
resulted from the same excess water being deposited on saturated ground. 
Observations (1953) on the impervious snow lysimeter at the snow laboratory 
in the Central Sierras near Soda Springs indicate that the permanent retention 
of liquid water in a 100-in.-thick snow pack of 31% initial density during a 
5-in. rainstorm was approximately 1 in., which included 0.3 in. as a result of 
the increase in temperature of the pack (isothermally) to 0°C. There was, in 
addition, transitory storage of liquid water in the pack, in the amount varying 
from 1 in. to 2 in., but this water drained out completely within 10 hr after 
cessation of rainfall and snow melt. In deriving the water-holding capacity of 
the snow, Mr. Riesbol assumed for Table 4 that the snow pack must reach a 
density of 40% before any water excess is available at the ground for contribu- 
tion to the stream flow. From the observations at the Central Sierra Labora- 
tory, runoff from the snow pack commenced after the pack density increased 
to 32.5%; the maximum density was 34%. In Table 4, it is not clear how the 
increments of snow melt shown in Col. 4 could total 11.0 in. when the initial 
water equivalent of the snow pack was 6.0 in. Also (considering the total 
storm period) there were 29.6 in. of rainfall and 6.0 in. of snow-water equivalent, 
of which 0.5 in. remained in storage in the pack at the end of the period and 1.5 
in. were retained initially in the ground. Thus, a total of 33.6 in. was available 
at the ground surface during the storm period. The total water excess, which 
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is the sum of the increments listed in Col. 13, is 26.63 in. The loss then is 7.0 
in. during the 54-hr perioil, which is equivalent to an average loss rite of 0.13 
in. per hr. For design-flood conditio1s, the value appears to be rather high 
for this area. In general, the retention and detention of liquid water in the 
relatively shallow snow pack shown by this example appears to be too great. 
Herzerr 8. Riespor,” M. ASCE.—The writer is indebted to Messrs. 


® Asst. Chf., Hydrology Branch, Project Planning Div., Bureau of Reclamation, U. 8. Dept. of the 
Interior, Denver, Colo. 


Gerdel, Hafterson, Linsley, and Rockwood for their constructive discussion and 
criticism of the paper. The purpose of the paper was to demonstrate the 
engineering application of scientific knowledge and data of snow hydrology to 
the planning, design, and operation of multiple-purpose reservoirs. Mr. 
Gerdel has urged that this application not include the use of hypotheses incom- 
patible with “* * * experimental evidence or the known physical properties 
of snow * * *.”” Mr. Hafterson emphasizes the fact that the application must 
be based on the judgment of the engineer, and that procedures must be con- 
tinuously developed to bridge the gap between physical knowledge and engineer- 
ing application. Each of the discussers had emphasized the need for new and 
improved methods of measurement and for more basic data. The writer agrees 
with these recommendations. 

Mr. Gerdel has shown that, under a variety of conditions, the density of a 
ripe, melting, actively discharging snow pack can vary over a wide range. 
Undoubtedly, this is correct. However, the writer has made the conservative 
engineering assumption (under the heading, ‘“‘Application to Spillway Design: 


Flood Produced by Rain on Snow’’) that ‘““* * * no release of water from the 


snow pack would occur until the snow attained a density of 40% 
This assumption was based on an examination of considerable data which has 
been corroborated by Walter T. Wilson, A.M. ASCE. Mr. Wilson has 
stated®; 


* “Storage and Movement of Liquid Water in a Snowpack,”” by Walter T. Wilson, U. 8. Weather 
Bureau, Washington, D. C., October 12, 1953. 


“The appended table gives average density as a function of number of days 
since establishment of the snow cover. The values given are based on a 
large mass of data including shallow and deep snow, short seasons and long 
seasons, high and low elevations and a wide range of other conditions. Of 
course, particular instances frequently vary considerably from the average 
values, Put experience has shown the average values to be reliable * * *.’’ 


Mr. Wilson’s table of the average-density trend contains the following 


information: 


Time since beginning of Average density 
snow cover, in days (percentage) 
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Most of the snow floods which are dealt with in the design and operation of 
multiple-purpose reservoirs occur within 100 days to 200 days after the begin- 
ning of the snow cover. 

Mr. Gerdel has stated that the temporary detention of liquid water in snow 
is a function of the time required for the development of internal flow channels, 
and Mr. Hafterson states that the detention varies with the depth and quality 
of the snow. In this connection Mr. Wilson® presents the following information 
and practical formula for water movement through a deep snow pack: 


“Tf a dense snowpack is deep enough to warrent routing, instead of merely 
lagging the rain or melt * * *, after its waterholding capacity has been 
satisfied the rates of rainfall or snowmelt may be expressed as a hydrograph 
of inflow to the pack at its upper surface. The outflow hydrograph at the 
snow-ground interface may be approximated by lagging the inflow and 
routing it by reservoir analogy. If the snow depth is d feet, the lag is 0.3d, 
and the routing coefficient is 0.5d, both in hours. Using the Muskingum 
notation, K in hours equals half the snow depth in feet, X is zero, and the 
hydrograph is translated a number of hours equal to three-tenths of the 
snow depth in feet.” 


Mr. Hafterson has shown that, in the derivation of a design flood produced 
by rain or snow, consideration should be given to the melting of snow by the 
heat of the rain water. As Mr. Hafterson states, this source of melting is 
usually regarded as being of minor significance ; however, in the Sly Park study, 
it appears that the melt from this source might contribute considerably to the 
flood runoff. Usually, the rain forms at a higher elevation and lower tempera- 
ture than is indicated by the surface dew point. Thus the rain would be sub- 
ject to some evaporative cooling during its descent, and, as a result, when the 
rain reaches the surface, it might be at a temperature several degrees lower than 
vhe surface dew point. This lower temperature would reduce the melt po- 
tential. 

Mr. Linsley notes the similarity between the hydrothermogram and his 
(1943) procedure and states that the method is basically the same as that 
offered by Mr. Clyde (1931). It is believed that any valid procedure for 
constructing the hydrograph of discharge from snow melt would necessarily 
resemble Mr. Linsley’s procedure. However, the only similarity between the 
writer’s procedure and Mr. Clyde’s procedure for net yield is in the assumption 
of a constant volume of melt per degree day. Actually, a more general and 
more flexible procedure is now (1953) in use in which both the basic temperature 
and the discharge factor are considered to be variable functions of the aceumu- 
lated percentage of the seasonal runoff volume. Methods have been developed 
by engineers of the USBR for constructing and adjusting the relationships 
among these functions, using only temperature and runoff data, so that a knowl- 
edge of the elevations, the area of the melting zone, and the lapse rate is not 
essential. Such data, if readily available, could perhaps be used to refine the 
procedure. Data concerning the wind velocity, humidity, and solar radiation 
could also be used to refine the procedure; however, these data are seldom 
available for design-flood studies. 

The hydrothermogram procedure developed by Mr. Barnes differs from the 
methods cited by Mr. Linsley as follows: 
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1. Ground-water flow and direct runoff are computed independently and 
concurrently from separate hydrothermograms. Different lag times, depletion 
factors, base temperatures, and discharge factors are used for the two elements 
of flow. 

2. The base temperature and discharge factor are considered to be variables 
and are determined from the accumulated percentage of the seasonal runoff. 

3. The runoff volume generated each day is expressed directly as the maxi- 
mum daily discharge in cubic feet per second. 

4. The subsequent distribution of the runoff is accomplished by combining 
all the flows (of direct runoff or of ground water) that are in the channel and by 
applying the appropriate constant depletion factor to the total. 


Mr. Rockwood is correct in stating that the methods presented by the 
writer for the derivation of the spillway design flood resulting from snow melt 
will not yield the maximum possible flood potential from that source. These 
methods are used as a rational procedure to obtain a maximum probable flood 
‘event which would have a magnitude somewhat less than the maximum possi- 
ble event. This maximum probable flood is of sufficient magnitude to justify 
its selection for prudent design purposes. Snow floods developed by these 
procedures are considerably in excess of maximum experienced floods for the 
area of derivation. 

Mr. Rockwood inquires how the increments of snow melt in Col. 4, Table 4, 
could total 11.0 in. when the initial water equivalent of the snow pack was 6.0 
in. The values in Col. 4, Table 4, are melt potential computed from Eq. 3 and 
are used only as a method of distributing the depth of snow on the ground at 
the beginning and end of each 6-hr time increment (Col. 5 and Col 6, Table 4) 
after accounting for the increase in density (Col. 7, Table 4) up to the 40% 
limit. Following the attainment of a pack density of 40%, the values shown 
in Col. 4, Table 4, reduce the waterholding capacity of the pack (Col. 8, Table 
4) by equivalent amounts. Mr. Rockwood states that ‘“* * * a total of 33.6 
in. was available at the ground surface during the storm period.” Actually, 
this amount should be 35.1 in. as the 1.5 in. shown at the head of Col. 12, Table 
4, is the retention which was assumed to have occurred during a previous event. 
This assumption provides a reasonable starting point on the retention-time 
curve. The average loss rate over a period of 54 hr of approximately 0.16 in. 
per br falls within the range of retention rates derived by the writer for the area 
from hydrograph analysis. 

The cooperative snow-investigation program, referred to by Mr. Rockwood, 
has contributed greatly to the physical theory underlying snow hydrology and 
to the practical applications of that theory. The USBR has maintained a 
close contact with that program, including the contribution of approximately 
four man-years of technical assistance. As the results of the program are fully 
analyzed and presented, they will become even more valuable to all engineers 
engaged in the planning, design, and operation of multiple-purpose reservoirs 
that receive their inflow from snow-fed streams. 

Corrections for Transactions.—Kq. 3 should be changed to read as follows: 


D = 0.00184 K U [(t — 32) 4 9.00578 (p — 6.11)] 
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